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Abstract 

Using tlie Baxter's T-Q relation derived from the transfer matrix technique, we consider 
the diagonalization problem of discrete quantum pendulum and discrete quantum sine-Gordon 
Hamiltonian from the algebraic geometry aspect. For a finite chain system of the size L, when 
the spectral curve degenerates into rational curves, we have reduced Baxter's T-Q relation into 
a polynomial equation; the connection of T-Q polynomial equation with the algebraic Bethe 
Ansatz is clearly established . In particular, for L = 4 it is the case of rational spectral curves 
for the discrete quantum pendulum and discrete sine-Gordon model. To these Baxter's T-Q 
polynomial equations, we have obtained the complete and explicit solutions with a detailed 
understanding of their quantitative and qualitative structure. In general the model possesses a 
spectral curve with a generic parameter. We have conducted certain qualitative study on the 
algebraic geometry of this high-genus Riemann surface incorporating Baxter's T-Q relation. 
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1 Introduction 



In the early seventies, R. Baxter proposed the method of Q-operator and the T-Q relation in 
his renowned solution of the eight-vertex model and the spin ^ XYZ chain in soluble statistical 
mechanics [Q, ^. Since then, the method has played a powerful mechanism up to nowadays in 
the two-dimensional exactly solvable lattice models, and the corresponding quantum spin-chain 
Hamiltonians. The quantum inverse scattering /algebraic Bethe Ansatz method developed by 
the Leningrad school in the early eighties [0, 0] systematized earlier results on two-dimensional 
integrable lattice models, and paved the way for the far-reaching effects in both mathematical 
and physical development in the past two decades. In the study of massless lattice sine-Gordon 
model, Izergin-Korepin and Tarasov |jl^ found, within the framework of quantum inverse scattering 
method, the ^-operator ( with C^-operators entries) which satisfies the Yang-Baxter equation for 
the six- vertex model's i?-matrix with the parameter q being the A^th root of unity: = 1, ( 
see the formula (|3|) of this paper). A slightly modified version of this ^-operator also appeared 
in the study of chiral Potts A^-state model Q. On the other hand, there are Hamiltonians of 
physical interest, exhibiting an intimate relationship with certain systems derived from the transfer 
matrix for a fixed finite size L while A^ varying. For L=3, the Hamiltonian, proposed first in 

and then investigated in a rigorously mathematical manner in can be interpreted as a 
generalization of the Hofstadter model fl], P], a renowned Bloch system with a constant external 
magnetic field. Through the quantum inverse scattering method, one can calculate the spectrum 
of the Hamiltonian by solving the (algebraic) Bethe ansatz equation. In |13| we formulate the 
problem from the algebraic geometry aspect, and consider it as the Baxter's T-Q relation on the 
spectral curve via the Baxter's vacuum state [|3|, |5| . In addition, a general scheme of diagonalizing 
the transfer matrix of a finite size L by means of the Baxter's T-Q relation ( or the Bethe equation) 
on the spectral curve was also discovered. At present, it is rather difficult to extract the explicit 
quantitative data for the spectrum problem by this approach, due to the complicated functional 
theory of the high-genus spectral curve. Nevertheless, when the spectral curve is totally degenerated 
into rational curves, we demonstrate in this paper that the polynomial equation derived from the 
Baxter's T-Q relation gives rise to the usual Bethe ansatz equation in physical literature. And for 
-L = 4, the transfer matrix produces the Hamiltonions of the discrete quantum pendulum and the 
discrete quantum sine-Gordon (SG) model, appeared in |6|, p^ . 

In this article we make a thorough study of the discrete quantum pendulum and SG model in 
the context of Baxter's T-Q relation by using the transfer matrix technique, ( for the Hamiltonians, 
see ( p!7[ ) ( pIsD of Sect. 2). Though the formulation via this approach can be made on Hamiltonian 
chains of an arbitrary finite size L, the mathematical treatment given in this paper for these two 
specific models takes advantage of certain special features presented only in the size L=4. The 
spectral curves, upon which the Baxter's T-Q relation is formulated, have a very high genus in 
general; indeed it is of the order . However, by examining their algebraic geometry properties, 
these spectral curves are found to form a family of algebraic curves covering the elliptic curves, 
so one could hope that the elliptic function theory would eventually play a role in the solutions 
of Baxter's T-Q relation. In the rational-spectral-curve case for a finite size L, the geometry of 
spectral curves becomes trivial. However, it still inevitably requires some subtle analysis of the 
Baxter vacuum state to extract the essential data on solutions of the polynomial equation, then 
to carry out the necessary algebraic study of a certain "over-determined" system of g-difference 
equations for a Ath root of unity g, a difficult problem for a general L. In the cases of discrete 
quantum pendulum and SG model, where the size L is equal to 4, we are able to obtain the 
explicit solutions of Baxter's T-Q polynomial equation by taking the special symmetric structure 
of polynomial functions into account. The results are complete from both the quantitative and 
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qualitative aspects, and the method provides a sound mathematical treatment on problems raised in 
P, 10 1 . The main advantage of our approach to these problems is to make use of the explicit Baxter's 



vacuum state as the precise form of "Bethe ansatz ground state", which was only postulated in the 



previous articles in literature, e.g. [10|. This approach to the eigenvalue problem of Hamiltonians 
has also been shown in ||l^ for L=3 on the Hofstadter-type models. It seems to us that the use of 
Baxter's T-Q relation would be more fundamental, and mathematically tractable, than the usual 
Bethe Ansatz technique. Furthermore, in this approach a certain mathematical theory of g-Sturm- 
Liouville type problem would possibly arise to entangle with the Baxter's T-Q polynomial equation 
originated from certain physical problems. 

This paper is organized as follows: In Sect. 2, we review the basic construction of Baxter's T-Q 
relation using the transfer matrix techniques, with the spectral curve of a high genus depending 



on the size L and parameters involved in the Hamiltonions. Some facts in |13| will be recalled 
here for the sake of completeness. Then we describe the constraints of parameters appeared in 
the spectral curves for the discrete quantum pendulum and SG model. In Sect. 3, we discuss 
a canonical procedure of reducing the Baxter's T-Q relation to a polynomial equation when the 
spectral curve is degenerated to rational curves for a finite size L. By converting the parameters to 



one special case studied in our earlier paper |13|, we obtain the Baxter's T-Q polynomial equation. 
In Sect. 4, we apply the conclusion of the previous section to the case L=4:, and with the further 
parameter constraints of the discrete quantum pendulum and SG model. The symmetric Baxter's 
T-Q polynomial relation is introduced, and a general discussion on the qualitative nature of its 
solutions is given. In Sect. 5, we construct explicitly the complete solutions of symmetric Baxter's 
T-Q polynomial equations, among which the discrete quantum pendulum and SG model for the 
rational spectral curve are included. Both the quantitative and qualitative nature of these solutions 
are revealed in the procedure of the derivation; in particular, the connection of the solutions with the 
usual Bethe ansatz equation in literature is clarified. In Sect. 6, we consider the discrete quantum 
pendulum and SG model for the general spectral curves, which are high-genus Riemann surfaces. 
By examining the geometric properties, we discover a canonical relationship of these curves with 
elliptic curves. A primitive qualitative analysis is made on the geometry of these spectral curves in 
connection with the eigenvalue problem of physical models through Baxter's T-Q relation. In Sect. 
7, we present the conclusion remark with a discussion of our future programs. We end with the 
appendix of presenting the identification of the sine-Gordon integral in |]l^] with the Hamiltonion 
given in this paper. 

Notations. To present our work, we prepare some notations. In this paper, Z, R, C will 
denote the ring of integers, real, complex numbers respectively, N = Z>o, Zjv = Z/A^Z, and 

i = V — 1. For a positive integers n, we denote C the tensor product of n-copies of the vector 
space C^. We use the notation of q-shifted factorials, 

(a;q)o = l, (a; q)„ = (1 - a)(l - aq) • • • (1 - aq""^), nGN. 



We shall minimize the repetition of materials in our previous article ]13| j, and so opt to use the 
same notation conventions as much as possible. 



2 Transfer Matrix, Baxter Vacuum State and T-Q Equation 

In this section we first recall some formulae in quantum inverse scattering method for later use. 
Most of the materials can be found in ||l^, including the original references. Then, we specify our 
discussion on the cases of the discrete quantum pendulum and SG Hamiltonian, the models we 
shall mainly consider in this work. 
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In this paper, N will always denote an odd positive integer with M = [y]: = 2M+1 (M > 1), 

and w is a primitive A^th root of unity, q := u>^ with = 1, i.e. q = co^^^^. A vector v in is 
represented by a sequence of coordinates, {vk \ k G Z), with the A^-periodic condition: Vk = Vk+N, 
i.e. V = {vk)j^^2ijsi' '^^^ standard basis of will be denoted hy \k) (A; G Z^v), with the dual basis 
{k\ {k e Zn) of C^*. 

Let Z,X he a pair of generators of the Weyl algebra with the A^th power identity : ZX = 
LoXZ, Z^ = X^ = /, and denote by Y{:= ZX) the element in the Weyl algebra. One has 
the relations, XY = uj~^YX, YZ = uj~^ZY, and Y^ = 1. The canonical representation of the 
Weyl algebra is the (unique) A^-dimensional irreducible representation, with the following matrix 
realization: 



Z{v)k = uj"-Vk , X{v)k = Vk^i 



Y(v) 



for V = (vk) G C 



N 



By using the operators X, Y, Z and the identity /, one can form a solution of the Yang-Baxter (YB) 
equation for a slightly modified i?-matrix of the six- vertex model, appeared first in [Q], and then 
studied in great details in |13]. The L-operator, depending on a parameter h = [a : b : c : d] in the 
projective 3-space P^, is given by the following 2x2 matrix with operator-valued entries acting on 
the quantum space C^: 

= ( xlz "^^f ) ' ^ ^ ^ ' 
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and it satisfies the following YB relation: 

R{x/x'){Lh{x) (g) 1)(1 (g) Lh{x')) = (1 (g) Lh{x')){Lh{x) (g) l)R{x/x') , (2) 

aux aux aux aux 

where the script letter " aux^^ indicates an operation taking on the auxiliary space C^, and R{x) is 
the matrix of 2-tensor of the auxiliary space with the following numerical expression, 



R{x) 



( 



xuj — x~^ 

uj{x — x^^) uj — 1 

LO — 1 X — X~ 





XUJ — X ^ j 



The operator (||) is related to the the following L-operator appeared in [|T2[ for the study of sine- 
Gordon lattice model by using the six- vertex model R-matrix and the Weyl operators U^V UV = 
q-^VU, = V^ = 1): 



/ xq — X 



V 







X — X 





q-q" 

X — X' 











(3) 



xq — X ^q ^ j 



with the YB relation: R(xlx'){{Ll{x)%\) ® (l{g)L^(x')) = ((10L^(a;') ® (L^(x) l))R(x/x'). 
Indeed, by the identification: Z = VU, X = V^^U, ( equivalently, U = q~Y2^V = q2ZY~,) 
Lh{x) in (|^) and L^{x) in (^) are related by 

LUx) = Lh{x)Y^q-2 . (4) 

By the matrix-product on auxiliary spaces and tensor-product of quantum spaces, the following 
L-operator of a finite size L, depending on the parameter h = {ho, . . . , /il-i) £ (P^)^, 

L-l 

^h(^) = (g (^) (— ^ho (x) L/,, (x) ... Lh^_^ (x)) , 

j=0 
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again satisfies the YB relation (|2[), hence its trace gives rise to a commuting family of transfer 
matrices: 

Tj^{x) = trauxiLj^{x)) , X £ C . 
The same conclusion also holds for Li(x), and T£{x), defined by 

L-l 

Lt{x) = (g) Ll^ix) , TSix) = tTau.{Ll{x)) . 

j=0 

By , one has the following relation of the above two types of transfer matrices : 

T;{x) = Tj^{x)D^ , where D := q-^ d Y. (5) 

For later use, we now summarize some basic facts on the Baxter's T-Q relation and the Baxter 
vacuum stateQ over the spectral curve in the diagonalization problem of transfer matrix Tj^{x) (for 
the details, see [jl3|). Applying the following gauge transform on LhAx): 



Lh^{x,^j,^j+i)=AjLh^{x)Ajl^, where ^j=( ] ^^'^^ ) (0<j<L-l), and := , 
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one has the expression: 



i.,^.ix,^„^,+ij - 1^ ^^^(^^^^. - 1) ) ■ 

Here the operator ^, is defined by F/j(x, ^, := F,'aY — xbX + ^'^xcZ — S,d. Accordingly, 

for h G (P'^)^ and = {F,o, ■ ■ ■ , £ ( C^)^, the modified L-operator becomes 

( Lr.AxJ) Lr,JxS)\ 
i^K(-.0:=(g)i..(x,e„e,+i)= '1 ,and ^ := ^0 • 

As the gauge-transform procedure keeps the trace unchanged, we have the relation: T^^ix) = 
tiCaux{Lj^{x,^)). For a given h, we will consider the variable (x,^) in the curve defined by 

cN N _ NuN 

We shall call the spectral curve in this paper. Over C^, we have the Baxter vacuum state defined 

L J. J 

by the following family of vectors in C : 

\p) := \pq) ... \pL-i) G i G C^) , 

where \pj) is the vector in governed by the relations: 

/0|„.\ = 1 (m|p,) _ ij+iaju:^ - xbj 



^The Baxter's T-Q relation and the Baxter vacuum state here were called by the Bethe equation, the Baxter 
vector respectively in |l^ 
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The constraint of the variables {x, ^) in ensures that the following properties hold for the Baxter 
vacuum state, 

where A±,t± are (rational) functions and automorphisms of C^, defined by 

L-l 

A_(x,^o, • • • ,^L-i) = Y[{dj - x^j+iCj) , 

A+(x,eo,...,eL-i) =yI ^^M-x%c,) ^ 

fJ-Q ij+iaj - xbj 

T± : (x, Co, • • • , Cl-^i) ^ {q^-^x, q'^^o, • • • , q~^^L-i) ■ 

This implies that, under the action of the transfer matrix, the Baxter vacuum state is the sum of 
its r±-translations: 

Tii{x)\p) = \t^p)/\-{p) + |r+p)A+(p) , for p G . (9) 

For a common eigenvector {(p\ C * of Tj^{x) (x G C), its eigenvalue A(x) is a polynomial of 
x, i.e., A(x) G C[x]; while its values on the Baxter's vacuum state, Q{p) := {^\p) {p G C^), is a 
rational function on C^. Then the following relation hold for f^{x) and Q{p): 

A{x)Q{p) = A_(p)Q(r_(p)) + A+(p)Q(r+(p)) , for p G , (10) 

which will be called the Baxter's T-Q relation for Tj^ix) Note that the {ip\ is also a common 
eigenvector of TS(x) with the eigenvalue A*(x) = g"A(x), where q'^ is the Z?" -eigenvalue of {ip\. 
The Baxter's T-Q relation for Tl{x) becomes 

A*{x)Q{p) = A*_{p)Q{r-{p)) + A+{p)*Q{T+{p)) , for p G , (11) 
with A*(x) = g"A(x) , A*_{p) = g"A_(p), and A+{p)* = g"A+(p) for n G Zat. For L=4, one has 
r^(x) = To + x^Ts + x^r4 , (x) = To* + x^Tl + x^T^ , (12) 

4 AT . 

where T2j are operators of (X" C with the expressions: 
To = aoaia2a3Uj'^ D + dodid2d3 , 

^2 = aoaib2CsY ^ X ^ Z + boCia2a3X (g) Z (g) F ® F + aobiC2a3Y ^ Z 
+aobid2C3Y X ® 1 Z + 60^10203^ 010 Z 0Y + bodid2C3X I 1 Z 
+{ajY ^ dj, bjX ^ CjZ) , 

T4 = boCib2C3DC~^ + cobiC2b3C , C:=Z®X0Z®X. 

One can also obtain the expressions of Tg^-s from (^). Note that the operators C,D T2 and 
commute each other. 

For the study of discrete quantum pendulum and SG model in this paper, we shall restrict 
ourselves on the case L=4 with the following constraint]^ in h, 

Qjdj = q^^, bjCj = —k^^ , for j = 0,2 , 

ttjdj = q~^, bjCj = —k , for j = 1,3 . (13) 



^The convention we use here is in tune with the one in |^. 
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where fc is a complex parameter. Then the operators T2j,T2j in (|12D now take the following forms: 
To = + dMd, , T, = T^DC-' + , 

+I&-,DU,~' + jSk^U^' + Ht^i + '-^DV,-' + T^V, + m W^-^ ; 
^^0= ,,^Dh^ + dod^d2dsD^^ , Tt=^^mC-^ + ^^D^C ,^ (14) 

' docirfs ^ kdodiC2 i ' did2C3 ^4 ' gd2C3 ' doci 1 

^kqdiC2^ + kdsco^ M ; 

where Uj,Vj are operators defined by 

C/i = Z0X0lOl, C/2 = 10Z®X®1, [/3 = 1®1®Z®X, C/4 = X®1(8)1(8)Z, 

yi = iozoyox, V2 = X(g)i(g)Z(g>Y, V3 = Y(g)X(g)i(g)Z, Vi = Z(g>Y(g)X(g)i . 

It is easy to see that the following relations hold among the above operators, 



Uj+iU, = uU.jUj+i , Vj+iV, = co^VjV.+i, (C/5 := Ui, V5 := Vi), 

UiUj = UjUi , ViVj = VjVi ifi=j (mod 2); 

Vi = U3U2 , V2 = U^Us, V3 = UiUa , Va = U2U1 ; 

U1U3 = C , U2U4 = C-^D , V1V3 = V2V4 = LoD . 

Under the constraint (p!3|), the functions A± in (^) become 
A-(x,^o, • • • ,^3) = nLo('^i - x^j+iCj) , 



(15) 



A c t ^ = dodid2d3{l+x'^qk''^f{l+x^qkf 

+ { ,?0,---,?3; (l+x^-idoc-'gfc-i)(l+x^-'d2C2'9fc-i)(l+^'Cj'(iic-'gfc)(l+x5-^ci3C3'gfc) 



By iQ), the Baxter vacuum state \p) = 0^=Q\pj) have the following expression, 

I cm , 

{m\pj) = 



^7-.^-(.g,;.c,i-^,^; for even J, 



In this paper we shall mainly study the diagonalization problem of , or equivalently T2 , under 
the condition (p^), plus the following further constraints on parameters dj,Cj and the operator C: 

dodM = 1 , = k^^'c^c^ , C = . (16) 

The above properties on the parameters arise from the connection of T| with the following physical 
models. 



(I) Discrete quantum pendulum. This is the situation under the constraint ([iq ), and with the 
further identifications: 

_ docithjj _ kcQ jj^i,, ^ . ci 1 _ kdQdiC2 . 

kc2 «ia2C3 kcod2d3 C3 



By (|15D, one has 



^ ^ doci -u. doci I I Cods coda i 

Qn-lQn = -, W ^Vx = , Qn-\Qn = — = — r^^2 • 

£12 C3 (22 C3 C2ai C2ai 
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Then in ( |14D become 

Tg* = r| = 2 

-T| =2{Q^ + Q~^ + Q^^i + Q-^_,) (17) 

+k{qQn-lQn + g^^Qn^Qn-l) + (^Qn^r^il + q'^Qn-lQn^) ■ 

The above —T2 is the Hamiltonian of discrete quantum pendulum in j^, subject to the following 
evolution equation: 

Qn+lQn-l = ( _|_ „'!3i )^ ' QnQn-1 = q'^Qn-lQr. 



(II) Discrete sine-Gordon (SG) Hamiltonian. This is the situation under the constraint (|T 
with one further identification: 



cidod3 ^ kco „i 

— D2U2 = — r^i?2[/ 

A;c2 C3dia2 



In this case, we have 



The above — can be identified with the discrete quantum sine-Gordon integral in |^], for which 
a detailed description will be given in the appendix of this paper. 



3 The Baxter's T-Q Polynomial Equation for Rational Degener- 
ated Spectral Curve 

In this section, we derive the Baxter's T-Q polynomial relation for a size L when the spectral 
curve is degenerated into rational curves, i.e., is a disjoint union of finite copies of the base 
x-curve. We shall reduce the general degenerated situation to one special case, which we have 
already discussed in our previous article [Ol. 



By the rational degenerated spectral curves, we mean the coordinates ^j^ of to be constants, 
i.e., independent of the variable x for all j. Then the parameters hj and the variables are subject 
to the relations: 

^3 "^J _ "i+l'^i+l f:2N _ for n < 1 < y 1 

^3 "-3 



In this situation, we define 



Ibj-idj-i 



i^T^L- (19) 



Then contains the following T±-invariant curve C, over which we shall formulate the Baxter's 
T-Q equation, 

C := {(x, eo, . . . , ih-x) I ^^0 = . . . = = ^^ ^ G Ziv} . 
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We shall make the identification of C with x Z^r via the following correspondence: 

C = P^ xZn , {x, rog^ . . . , rL_iq^) < > {x, I) . 

Then the automorphisms t± on C are expressed by 

r± : {x,l) ^ {q^^x,l- 1) , 
and the action T(x) (:= T^{x)) on \x,l) in (^) now takes the form: 

T{x)\x, I) = \q-^x, I - l)A_(x, I) + \qx, I - l)A+{x, I) , (20) 
where A± are the following rational functions of x: 

L-l 

A_(x,/) = {do---dL-i)Yl{l- xq^dJ^Cjrj+i) , 

j=0 

1 - x'^a-^d'^h.c^ 

A+(x,o ={do---dL-i)X{- ■ 

With the substitutions, 

(4 • • • (iL-i)"^T(x) ^ r(a;) , (do • • • dL_i)-iA±(x, /) ^ A±(x, , 
the relation ( pO[ ) still holds for the modified A±, now with the expressions: 

L-l L-l 1 _ ^2^*2 



A_(x,0 = \{{l-xc]q') , A+(x,0 = n T 

i=0 j=0 J ' 

where 

c* := dj'c.rj^, (= a7^^rT^\) • (21) 

Furthermore, one can convert the expression (|7|) of the Baxter vacuum state over C to the following 
component-expression of the Baxter's vector \x,l): 

(k|x,/) = ^1 I [[ 



{^c*q^+^;u;)k, 

Here the bold letter k denotes a multi-index vector k = (/cq, . . . , A;l_i) for kj G Zjv, and the 
square-length of k is defined by |kp := ^2^=0 k'j. Each ratio-term in the above right hand side is 
given by a non- negative representative for each element in Zat appeared in the formula. With the 
above description of T{x) on the Baxter vacuum state \x,l), the discussions of Sect. 4 and Sect. 
5 Proposition 2, 3 in can be applied to our present situation. This enables us to state the 
following result on the Baxter's T-Q equation and its connection with the transfer matrix T{x): 

Theorem 1 Let f^,f° be functions on C defined by 



r(x,2n)= n ^ r(x,2n + l)=n ' 



n+l 

ixc*;u;)n+i ' ' ' ^' M, {xc*q;uj)n+i 
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For X G P and / € Z^r, we define the following vectors in ® C , 

N-l 

|x)f = ^ \x, 2n)r (x, 2n)J^ , \x)f = En=o k> 2n + 1)/° (x, 2n + 1 V" , 

n=0 

k)/" = k)f9~''"(92;) + \x)fu{x) where ■u(x) := Hj'Jo^ll - c*^){xc*q;q'^)M ■ 

Then 

(i) |x)fn(q'x) = \x)°q^u{x), or equivalently, \x)'l~ = 2q~^\x)'(u{qx) = 2|x)°ti(x). 

(ii) The T{x) -transform on \x)l' is given by 

g-'r(x)|x)+ = |(7-ix)+A_(x,-l) + |(7x)+A+(x,0) , I £ Zn ■ 

(Hi) For a common eigenvector {ip\ of T{x) with the eigenvalue A(x), the function Q'j^{x){:= 
((^|x)^) and A(x) are polynomials with the properties: 

degQ+(x) < (3M + 1)L, degA(x) < 2[-§], A(x) = A(-x), A(0) = g^' + i, 

and they satisfy the following Baxter's T-Q equation: 

L-l L-l 

g-'A(x)Q+(x) = n (1 - xc*q-')Qlixq~^) + J] (1 + ^c*j)Qti^<l) ■ (22) 
j=0 j=0 

Furthermore, for < m < M , Q^{x) and Q^_^{x) are elements in x™ 11^=0^(1 ~ x^ c*^) C[x]. 

□ 

For the rest of this paper, the letter m will always denote an integer between and M, 

< m < M . 

By {iii) of the Theorem the equation ( p2| ) for the sectors m,N — m can be combined into a single 
one by introducing the polynomials Am(x), Q{x) via the relation, 

L-l 

(A^(x), X- n (1 - x''cf)Qix)) = {q-^K{x), Q+ (x)), ((/"A(x), Q+_„(x)) . 

Then the equations ( p2| ) for / = m,N — m are equivalent to the following polynomial equation of 
Q(x), Am(x): 

L-l L-l 

A„(x)Q(x) = n (1 - ^c*g-i)Q(xQ-i) + H + ^c*)Q(x(?) , (23) 
i=o i=o 

with the following constraints of Q{x) and Am(x), 

deg Q(x) < ML - m, deg A^(x) < 2[^^], A^(x) = A^(-x), A„,(0) = + (7"'". 

By (|), the above A m(x) is indeed the eigenvalue of T*(x); while ( p3[ ) corresponds the Baxter's 
T-Q equation (|Tl| ) for T*(x) on the sectors m, N — m. 
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For L even, by the construction of Ti{x), one can see that the eigenvalues of T£ are non- 
zero, hence deg Am{x) = L. In the situations we will consider later on, the polynomial solutions, 
Q{x), Am{x) of the equation ( p3D possesses certain reciprocal symmetry property. Here we call a 
polynomial P{x) to be reciprocal if P^{x) = P{x), where P^(x) is the polynomial defined by 



Proposition 1 For L even, assume that the polynomial h.m{x) in (|22| ) is reciprocal, and the pa- 
rameters CjS and the degree d of Q[x) satisfy the following properties: 
(i) {cl, = {-c^-^g, . . . , -ct_\q} . 

Then Q^x) is also a solution of ( p^ for Am{x). 

Proof. By substituting x by x^^ in (22), and then multiplying x"^"*"^ to the equation, one obtains 
the relation, 

Am{x)QHx) = 9™+' n n (1 + xc*-')QHxq-') + '^-^ n nV - xc*f'q)QHxq) , 
j=0 j=0 j=0 j=0 

By (m), we have 

L-l L-l 
q"'+'^ n = Q'"" , q-"^-d~L -Q ^* ^ 
j=0 j=0 

Then, by (i), the above equation of Q^{x) is the same as (23). □ 



The following algebraic fact was shown in |13] Lemma 6, which we just state here for later use. 



Lemma 1 Let n he an odd positive integer, A be anx n-matrix with complex entries aij satisfying 
the relations 



a. 



[-iy+^+'^an-j+i,n-i+i , for 1 < i,j < n . 



Then A is a degenerated matrix. 

□ 

4 The Baxter's T-Q Polynomial Relation for L=4 

For L=4 rational degenerated case, the parameter h we discuss later in this paper will subject to 
the constraint ([l3|), and be confined only to the following situation: 

J 1 L f for even j, 

qaj =dj = l , -bj =cj = \ . (24) 

[ /c2 tor odd J. 

Then, by (|l9|) and ( ^T|) we have rj = (—q)^ for all j , and 

c* = I (-9)^^^ for even j, ^25) 
1 (— g)2A;2 for odd j. 
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The operators ( [L4[) in Ti{x) become 



t: = d^c-^ + d^c 



-1 -1 -1 -1 , -1 , -, -1 



= D—Ui + D—U2 + D—U3 + D—U4 + kq-^D—Vi + k-^q-^D—V4 



By = 1, the polynomial Am(x) in (p^) now takes the form: 



0<m< M, 0</<2M 



(26) 



where A is an eigenvalue of — r| for the sectors {m,l), {N — m,l), which can be regarded as the 
label of the eigenvalues of D~ and C^^. The Baxter's T-Q polynomial equation (23) now takes 
the form, 



with 



'A{xi-q)-)Qixq-') + q"'A{x{-q)2)Qixq) 

1 1 



(27) 



A(x) = (1 + 2cx + x^)^ , c:=^{k^+k^) 



and deg Q{x) < AM — m. The polynomials Am,i{x),Q{x) will be called the eigenvalue and the 
eigen-polynomial of (^) respectively, whenever Q{x) is a non-trivial function. Indeed, the value A 
in the expression of A^^i^x) is an eigenvalue of — T2 . In the following, we are going to study the 
equation (^) for a generic A in (p6|). 

For the rest of this paper, the parameter c will always be a generic complex number unless 
otherwise stated. The polynomial Q{x) will be denoted by 



Q{x) = J2^ 

j=0 



and we define a 



d := deg Q{x) , 

J .— for j not between and d. An equivalent formulation of the equation (|27| ) 
tem of difference equations in A and ajs : 

(28) 



is the following system of difference equations in A and ajs : 

ujaj + VjUj-i + {5j - A)aj_2 + ujajs + fijaj-i = , (j > 1) 
where the coefficients in the above equations are defined by 



jm+j _|_ q-m-j _ qin _ g-m ^ 

-(4c2 + 2)(g"^+J-i + 



-4ci(g 



m+j- 



4ci(g"+^-2 - 

: qm+j~2 _j_ -m-j+2 



(29) 



m+l 



-m—l 



Indeed, for the system (|28|), it suffices to consider those relations for the index j between 1 and 
d + 3. Note that the relations in (pSD for 2 < j < d+2 give rise to the following eigenvalue problem: 





/ ^d+2 


Ud+2 


fJ-d+2 





















Vd+1 


^d+l 


Ud+1 


fJ-d+l 




















Vd 


Sd 


Ud 


fJ-d 
































- A 




















> 



































Sa 
























S3 


U3 






V 




I 









V2 


V2 


S2 


J 


> 





ad-i 



(30) 



V "0 / 
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Hence for a solution of the system (^), the A can be regarded as an algebraic function of c, and it 
has a limit as c tends to some special value cq. 



Lemma 2 For equation (|27| ) with a given a (no generic property required), the degree d of Q{x) 
satisfies the following conditions, 

and the zero-multiplicity of Q{x) at the origin is equal to one of 0,N,N — 2m, 2N — 2m. 



Proof. The upper bound 4M — m of d is given by the assumption of (|27|). If d = 0, then a non-zero 
constant is a solution Q{x) of (p7|), and we have 



AmA^) = g-'"A(x(-g)^) + g-A(x(-g)2) . 

By the even- function property of Am /(x), the above relation implies g^m+i _ g2m+3 _ hence 
q'^ = 1, a contradiction to the odd assumption on the integer N. Therefore d > 1. Comparing the 
coefficients of the highest degree of x in (p7|), one has 



This implies q^+'^+'^ = q^+i qj- gi^'"^'^ i.e., ^"^+2 = qi or g-2m-'_ Denote r the zero- multiplicity of 
Q{x) at X = 0. By comparing the coefficients of degree r in (p7|), we have 

hence r = 0, —2m (mod N). Then the conclusion of r follows from d < 4M — m. □ 

Remark. For the results obtained later in this paper on certain special cases, and also on the 



similar problem of size L = 3 in |13], the solution Q{x) in ( |23D always possesses the property 



Q(0) 7^ 0. In this situation, one can write 



d ^ 

Q{X) = \[{X--) , Zj^Q. 

Substituting x = j- in (p7[), one obtains the following relations of ZjS : 
which is the Bethe ansatz equation appeared in literature, e.g. [Q. □ 

A special case in the above setting happens when Kj^Ax) in (p7|) is a reciprocal polynomial. For 



the convenience, the Baxter's T-Q relation (27) will be called a symmetric T-Q polynomial relation 
if the following condition holds: 

^m,«(^) = equivalently q^ = l,q~'^'^, i.e., l = 0,N-2m (mod iV) . 

In this situation, (|2^ ) becomes 

Am,l{^) = g"^ + q-"^ + + (g- + q-^)x' , (32) 
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and the coefficients (l29| ) in the system ( psD have the following form, 

^. ^ qm+j q-m-j _qm_ ^-m ^ ^. ^ 4ci(q™+i-| _ g-'^-J + l), 

Sj = -(4c2 + 2)(g'"+J-i + (33) 



3 ™ „■ I 3 



U 



Note that by the equalities, Uj+i = —vj and = I'j, the transport of the square matrix in pO| ) 
is unchanged after substituting c by — c. Hence the eigenvalue A for the symmetric T-Q polynomial 
relation necessarily becomes an algebraic function of c^, equivalently, the following property holds 
for A: 

A = A(c) = A(-c) . (34) 

Furthermore, the relation (p7| ) is unchanged when substituting (c, x) by (— c, — x); this implies that 
if Q{x; c) is a solution of ([27|) , so is Q{—x, ; — c). 

We now determine the qualitative nature of a solution Q{x) for the symmetric polynomial T-Q 
equation. 



Lemma 3 For a symmetric T-Q polynomial relation (27), there is no non-trivial solution Q{x) of 
degree d = N — 2 with the zero-multiplicity at x = equal to N — 2m. 

Proof. Otherwise, one has m > 1 and 

g(x) = j;^-2'^g(x) , where Q(0) / , degQ = 2m - 2 . (35) 
Write Q{x) = Y^'^=q'^ ajx^ . Then Q{x) satisfies the relation 

A^m,l{^)Q{x) = q^A{xi-q)^)Qixq-') + 9-"A(x(-g)5)g(xg) , (36) 

or equivalently, the coefficients ajs of Q{x) satisfy the following system of equations, 

i^jOj + VjOtj-i + {5j — X)aj-2 + Ujaj^3 + fljOtj^^, = , 1 < j < 2m + 1 , (37) 

where i'j,Vj,6j,Uj,Jlj are expressed by the similar forms as in (^) by changing m to — m in the 
corresponding term. By Vi ^ 0, we have m > 2. By the equalities, 

Z^j = ^2m+2-ji = U2m+2-j, = '52m+2-j j 

Q''^(x) also satisfies the equation (|36|). In general, for a polynomial (^(x) of degree d satisfies (|36|), 
(i = 2m — 2, N — 2 (mod A^), and the minimal possible d is 2m — 2. Hence the dimension of 
the solution space of Q{x) with degree < 2m — 2 is equal to one. Fbr Q(x) in (|5|), Q'f^x) is 
a scale- multiple of Q{x), which implies Q\x) = ±Q{x). Therefore, Q{x) is determined by the 
coefficients aj for < j < m — 1, which involve only those equations in ( |37|) with 1 < j < m + 1, 
subject to one of the following two conditions: aj = a2m~2~j for all j, or aj = —a2m-2~j ^or 
all j. On the other hand, the relation (|3^) is the same when we substitute {c,x) by (— c, — x), 
hence Q{—x, — c) = Q{x,c). Hence we may assume the coefficients cij = aj{c) satisfy the following 
properties: 

Qo(c) = 1 , Qj(— c) = {—iyaj{c), for all j . (38) 

For a solution {A, cij} of ( |37| ) for a generic c, A is a solution of the eigenvalues problem arisen from 
those relations for 2 < j < m. Hence A = A(c), an algebraic function of c so that the limit of 
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denoted by Aoo, exists as c ^ oo. For 1 < j < m — 1, by z^j 7^ one can conclude aj{c) = 0{c^) as 
c — > 00. Denote 



limM£), 

c— >oo cJ 



< J < m — 1 , 



and 



m+fc— ^ „nn—k+l 



for 1 < A; < m+1. By multiplying c on (^), and then taking the c-infinity limit for 1 < j < m+1, 
one obtains the following matrix relation on Aoo and ajs, 



/a' — \ 









J' — \ 

"m ^00 



V. 



m—l 

















V 





























% 












% 














( am-1 \ 

am-2 



. (39) 



ao 



Note that ao = 1, and 6'j 7^ 5'^ for 2<j^k<m + \. The square matrix by deleting the last row 
in (|39[) becomes the eigenvalue problem with m distinct eigenvalues, hence 



An 



5'i , for some 2 < / < m + 1. 



This implies = for / — 2 < /c < m — 1, and ai-2 / 0. In (|39[) , the row containing — Aoo 
gives the following relation of a/ and a;_i: 



(40) 



One would expect the right lower square matrix of size / — 1 with Aoo = S'l has non-zero determinant. 
This statement is valid for a small number / by direct computation, hence it leads to a contradiction. 
However, it is a difficult task to obtain a mathematical proof of such a statement for a general I. 
For our purpose, we are going to provide another way to justify the conclusion of the lemma using 

m + 1, then a^-i / 0. This implies Q^{x) = Q{x). The 

(41) 



(37). First we consider the case I 

relation (Hol) and the (m + l)th equation in (|37|) become 



i^m+l - A(c))am_i + 2Um+iam,-2 + ^Jj-m+ia.m-3 



. 



By the property (^) on A, the c-infinity limit of 0"^™"^-* -multiple of the second relation in ( |4l| ) 
gives rise to the equality, 

— 4am-i + 2u^_,_]^am,-2 = , 

which is incompatible with the first relation of (p]). Hence, we may assume 2 < I < m. By (38) 
and afc = for I — 2 < k, one has 



a'^jC*^ ^ -|- lower order term , as c ^ oo 



l-l<k<m-l. 
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By Aoo = ^[i one can take the c-infinity limit of c ^' -multiple of the lih. relation in (| 
yields the following identities: 



which 



when I < m - 
, when I = m. 



1; 



(42) 



Similarly the c-infinity limit of c ^^-multiple of the (/ + l)relation gives rise to the following 
relations: 



^i+i^'i+i + '"i+i^'i + i^i+i - W^i + = , when I < m - 2; 

i + i^'m- ^'m~i)^'m-2 + u'^^m-s = , when I = m - 1; 

when / = m. 



v' a' 



(43) 



m+l 



61 



m/°m—l 







For / = m, by using — ita772— 2i ( |4Q| ) ( [43[ ) and the last relation in (|42| ) will lead to a contradiction. 
For / < m — 2, we continue the same procedure on the c-infinity limit of c^*-multiple of the (s + 2)th 
relation for s > I, then obtain 

^'^+23^+2 + ^''3+2^3+1 + i6's+2 - KWs = • 

Hence one has the following relations for a^s, 



m+l 

v' 




5', 





6L - 6, 



rr' 



^m— 1 






A' 







6', 








^1+3 



''1+3 



V 



'1+3 

v' 



S', 



( < 
a 



: ■• ■• vi+2 ^1+2 "1+2 

which implies a^ = for s > I. The relations (42) (43) become 



/ 

m— 2 



^(3/-! + -^^1-2 + u'i3i^3 = , {S'l^i - d'i)a'i_i + nJ+ia/_2 



together with (40), this provides a contradiction to ai^2 7^ 0. □ 



For a symmetric ( p7| ) equation, the eigen-polynomial Q{x) has the following property: 

Theorem 2 Assume there exists a non-trivial polynomial solution for a symmetric T-Q polynomial 
relation ( p7|) with a given reciprocal polynomial Am,i{x). Then the equation has one- dimensional 
solution space, generated by a monic polynomial Q{x) of degree 2N — 2 — 2m with Q{0) 7^ and 
Qt(x) = ±Q{x). 

Proof. Let Q{x) be a non-trivial polynomial solution with the degree d. By Lemma |2[ we have 
d = N — 2, N — 2m — 2, 2N — 2m — 2. First we are going to show that d = 2N — 2m — 2. Otherwise, 
d is one of two odd integers, N — 2m — 2 ov N — 2. By Lemmas ^ and ^, we may assume Q(0) 7^ 
with ao = 1. The coefficients q^s of Q{x) satisfy the relation (p8|). When d = N — 2m — 2, we have 
Vj 7^ for 1 < j < d. This implies the x-coefficients Uj of Q(x)(= Q{x]c)) are polynomials of c 
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and A = A(c), hence aj = aj{c). As c tends to 0, the coefficients Oj = aj{0) of Q{x;0) satisfy the 
corresponding relation (p8|): 

UjOj + {6j — X)aj-2 + HjCtj-A = ,l<i<d + 3, 

with ai = 0. Hence = for odd j, and the polynomial Q{x; 0) has an even degree < N — 2m — 2, 
which impossible by Lemma It remains the case when d = N — 2 with Q(0;c) ^ 0. Now the 
dimension of the (5(x)-solution space of ( p7| ) with deg Q{x) < iV — 2 is equal to one. As Q{—x, ; — c) 
is also a solution of the symmetric T-Q relation, we have Q{—x,;—c) = Q{x;c), equivalently , 
aj{—c) = {—ly aj{c) for all j. Therefore Q(x;0) is again a polynomial in x with an even degree 
< N — 2, which contradicts to Lemma ^ for c = 0. Hence we have shown that any solution Q{x) 
with the eigenvalue Am,i{x) must have the degree d = 2N — 2 — 2m, which implies the (5(x)-solution 
space is of dimension one. By (p5[), the polynomials, Am,i{x) and Q{x), satisfy the conditions of 
Proposition H, hence Q^{x) is a also a solution of (27). Therefore Q'^ix) = ■yQ{x) for some non-zero 



constant 7, which implies 7 = 1 and Q{0) 7^ 0. Then the conclusion follows immediately. □ 

Remark. For a polynomial Q{x) in the above proposition, the roots of Q{x) are all non-zero; 
furthermore if x^ is a root, so is x^^. Hence the collection of all roots Xk (counting multiplicity) is 
the same as that of x^^s. The criterion for Q^(x) = —Q{x) holds if and only if Q{x) has the root 
X = 1 with a positive odd multiplicity. □ 



5 Solutions of Discrete Quantum Pendulum and Sine-Gordon Model 
in the Rational Degenerated case 

In this section we are going to derive the complete solution of symmetric T-Q polynomial relation 
(27); hence we now only consider the sectors (m, /) = (m, 0), (m, N — 2m). By Theorem ^, we may 
assume 

d = 2N -2-2m , Q^x) = ±Q{x). 



Hence the coefficients in (|33| ) possesses the following symmetric relations: 

The system (^) is equivalent to the eigenvalue problem (|30| ) together with one more constraint: 

uiai + viao = , (44) 

and the a^s satisfy either one of the following conditions : 

ai = Ud-i for < i < d, i.e., Q\x) = Q{x); (45) 
ai = —ad-i for < i < d, i.e., Q'^{x) = —Q{x). (46) 

Note that the polynomial Q{x) is determined only by the first half of its coefficients, i.e., ao, . . . , . 

2 

In the case (Bq), one has ad = 0. Furthermore, through the transformations: 

2 

Vj , Vj , 6j , Uj , fij fj^ji , Uj' , 6j' , Vji , Uji , where j' := d + 4 — j ; 
ai ad-i {0 < i < d) or ai —ad-i (0 < i < d) , 
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the equations for j > | + 3 in 



follow from those for j < ^ + 2. So we need only to consider 



the relations for 1 < j < ^ + 2 in (p8|), which are regarded as equations of A and for < A; < ^. 
Note that the + 2)th equation in ( pSD has the form : 



+ ad 

2 ^ 



= 



(47) 



which is a trivial relation in the case (|4q). 

By (p^), the rational degenerated case of discrete quantum pendulum and discrete sine-Gordon 
corresponds to C = 1, i.e. the sectors with / = in the symmetric (27) relation; in particular, by 
([TtI) and (|2^), the discrete quantum pendulum is given by D = C = 1, i.e. (m, /) = (0, 0) . 

Theorem 3 For the symmetric T-Q polynomial equation (^), there are N distinct eigenvalues X, 
each of which has one- dimensional eigenspace generated by a monic eigen-polynomial Q{x) of degree 
d = 4M— 2m with Q{0) andQ\x) = ■iiQ{x). Furthermore, there are (M+l) eig en-polynomials 
Q{x) of the type Q'^{x) = Q{x), and the rest M ones are of the type Q^{x) = —Q{x). In particular, 
the Baxter's T-Q polynomial relation of SG model are those for the sectors {m,l) = {m,0), and the 
discrete quantum pendulum is the one for {m,l) = (0,0). 



Proof. The relation (^41) is a non-trivial constraint for 0<m<M — Ibyz^iT^O; while for m = M, 



both ui and vi are zeros, hence (44) is a redundant one. In this proof, we shall first consider the 
case with m = 0, then 1 < m < M — 1, and finally on m = M. 

(I) m = 0, i.e. {m,l) = (0,0), which is the rational degenerated case of discrete quantum 
pendulum. We have d = 4M. Consider the relations for 1 < j < ^ + 2 in the system (28) as 
equations of A, ao, ad. By i^d,-, = 0, the problem is formulated in the following matrix form: 





Vd^ 


Sd,, — X 


Ud, 
2^ 


-1 


2^ 




Ud 

2 


Vd 
2 


5d - 

2 


A 


Ud 
2 







Ud 1 
2 ^ 


Vd 
2 


1 


6d -I - 

2 ^ 
















V 



together with the constraint 



iq + q~ 



2){ad_ 

2 







fid 

2 



\ 



1/4 












Vi 


64 -X 










63 -X 











82 -X 









Vl ) 



which is now given by 



( ad \ 



ad 
2 ' 



(48) 



+ 0^-2) + 4ci{q2 — q 2 ){ad,^ -\- ad 



V "0 / 



ic^ + 4 + X)ad = 



(49) 



Note that the square matrix of size (= | -|- 1) in ( [48D satisfies the condition of Lemma ||. Hence, 
by I'j ^ for 1 < j < ^, the system 
generated by a basis element (a^) 



0<fc<i 



has the one-dimensional eigenspace for any given c and A, 
with ao = 1- In fact, for 1 < /c < |, can be expressed 
by a polynomial of A and c, regarded as a polynomial in A with coefficients in C[c], and denoted 
by C(k = Pk(,X). The A-degree of Ofc is given by deg Pk{X) = [|]. In the case (p5D, the relation (HT 
becomes: 



2(g + q'^ - 2)pd_^{X) + 8ci{q2 - q—)pd_-^{X) - {8c^ + 4 + X)pd{X) = , 
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which defines A as an algebraic function of c. As the A-degree of the above relation is equal to 
M + 1, there are (M + 1) A-values for a generic c. In the case (^6|), ( [49D is a trivial relation. The 
relation ( ^8|) becomes: 



/ '5d_^-^ — A ^d+i 



2 



I'd 



5d — \ 

2 



2 ^ 

Ud 
2 












1 (5. 



A Uc 



^1 



V 



1^4 



f4 






f3 





li4 

53-A 





/i4 

<52-A 
Vi 



( "^-1 \ 



. 



Solutions of the above equation can be obtained from the system (^) alone, and then imposing 
the following constraint on A, 

ad = pd{X) = . 



As the A-degree of is equal to M, the above equation gives rise to M eigenvalues of A, with the 
corresponding eigenvector having the components Uk = PkW, 1 ^ k < f — 1- By Theorem ^, the 
Q(x)-eigenspaces are all one-dimensional, hence the conclusion follows immediately. 

(II) 1 < m < M - 1. We have I = 0, N - 2m, and d = 2N - 2 - 2m. By (||), for 1 < j < f + 1 
one has 



44> j = n := N - 2m 



(50) 



In the case (^5|), the relations for 1 < j < 2 + 2 in the system (28), considered as equations of A 
and ctfc (0 < A; < ^), can be formulated in the following form, 



S T 
u 



ai 

V ao y 



/ ad \ 



V On / 



(51) 



and together with the constraint (|47|). Here S,U are square matrices of the size (| — n + 1), n 
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respectively, T is the — n + 1) x n matrix, with the following expressions: 



+1 



2 ' 



2 



Vd 
2 



6d 

2 



A Ud 

2 



fid 

2 









V 

Vn-1 Sn-1 - A Un-l jJ-n-l 





^n+3 ^n+3 "n+3 

— A Un+3 

'5n+2 " A 

Vn+1 J 
\ 

: 



V 
/ 







^3 



V3 S3 - X U3 



/^n+3 
Un+2 
V Sn+1 - A Un+1 Hn+1 



U2 

••• 0\ 



V2 52 - X 
Vl Vl ) 




/ 



Note that there is the term Vd in the second entry of the first row of 5; while /Un+2 = in the 

2 

matrix T. By (p3|), the matrix [/ satisfies the condition of Lemma [l| for any A. Hence, by (|50|), 
the system = inside the relation (51) has the one-dimensional solution space generated by a 
vector with 

k 

ao = 1, Ofc =Pfc(A) G C[c][A], where deg pk{X) = [-] for k < n. 
Using the above vector -0, we consider the system 

sip = -Tip . 



By (50), one can first solve {k > n) in terms of A, c with the form : 

ttfc = rk{X)an + qk{X), where rk{X),qk{X) e C[c][X] , deg rk{X) + = deg qk{X) = [|]. 
Furthermore, satisfies the following relation : 

= Vd^-^Od + ii^d_^_-^ + Sd_^_^ — A)Qd_-^ + Ud_^_-^ad_2 + ^d_,_j^ad_3 = r(A)an + g(A) , 



20 



where r(A),g(A) G C[c][A] with deg r(A) + ^ 



0!k 



r(A) 



- deg q{X) = [!^]. Hence 
, Pk{X) := r{X)qk{X) - rk{X)q{\) . 



By multiphng the ak by r(A), we obtain a solution of for all A with the new afcj < A; < ^, in 
the form : 

a, = P,{X), deg P,(A) = [^] + [^] - ^ ; (52) 
in particular, deg Pd(X) = M. Now the constraint ( p7|) becomes 

2 

by which one can show that the above relation gives rise to (M + 1) A-values for a generic c. 

In the case (^6[), ( ^) becomes a trivial relation. We consider the following eigenvalue problem, 
similar to the one in ( pll) by changing 5 to S~: 



S- T 
U 



, 



ai 



/ Oid \ 



V an / 



(53) 



where the matrix S differs from S only on the (l,2)-th entry by changing I'd,, to — z^d , , , i.e.. 



5- 



+ 1 



2 



2 



(5d — A Ud 
2 2 





V 



Un+S Vn+3 (5n+3 " A Un+3 

l/n+2 '^n+2 " A 

I^n+l fn+1 



Then the problem (p^) with the condition ad = is equivalent to the problem in the case (|4q). 
As in the discussion of the eigenvalue problem (|5l|), there is a solution of (51), a^s, with the form 
Ok = Pk{X) satisfying the property (^). Then PJ is a degree M polynomial of A, and its zeros. 



"d = Pd (A) = 0, give rise to M A-values for the case (46). The conclusion now follows from 

2 2 
Theorem 0. 



(Ill) m = M. We have / = 0, 1, and d = N — 1. As the values of z^i, ui, un+2, in (29) are 

all zeros in this case, the relations of j = l,(i + 3 in (28) are redundant. Hence the system (|28| ) 
is equivalent to the eigenvalue problem (^0[) with d = — 1. In the case (|45|), the collection of 
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CKfc, < A; < |, among the coefficients of Q{x) is the solution of the following eigenvalue problem : 


















\ 


> 


























l^M VM 


UM 






























-A 
















> 























: ■ ■ . 1/4 




Sa 




















S3 








V 




I 





1^2 


V2 


<52 


) 







. 



V "0 / 



Note that the coefficients in the first and second rows have some extra terms, comparing to the rest 
of entries. There are (M + 1) A-eigenvalues for the above relation, which gives rise to the solutions 
for the case (H). The rest M A- values for (H) are those for the case (Hq). Then the conclusion 
follows from Theorem 0. □ 



Remark By Theorem ^, an eigen-polynomial Q{x) of the symmetric T-Q polynomial equation 
(]27|) are in sectors (m, /) = (m, 0), (m, N — 2m), and it has the form. 



1 



2Af-2-2m 

Qix) = n (x--), z,^0 

j=l ^3 



The Bethe ansatz equation ( |3l| ) of z^s becomes 

' ^ n — ' 1 < j < 2iV - 2 - 2m . 



By the description of Q{x) in Theorem |3|, one has the reciprocal constraint on Zj's, i.e., {zj}'^^-^ ^ ^"^ 
{z^^l^^f ^"^"^ (counting the multiplicity). □ 



6 The General Spectral Curve for Discrete Quantum Pendulum 
and Discrete Sine-Gordon Model 

In this section, we are going to explore the geometrical structure of the spectral curve (P) for 
the discrete quantum pendulum and SG model. By ( |l3|) and (|l^), the parameter h for the curve 
has the following constraints, 



^J.N ^ - - — (54) 



aj = q'^dj\ bj = -k-'^c-\ (ej := (-1)^), < j < 3 
dodM = 1, cf = e^c^c^ 



By the discussion in Sect. 8 of |jl^, the value ^j^s of the curve are determined by x^ , which 
we denote by y := x^,rj := ^j^. The variables {y,r]) satisfies the following equation of the curve 
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where Aj^, Bj^, C^, are polynomials of y given by the relation, 

In fact, by computation, one has the following expressions of these polynomials, 

A^^{y) = -(5-1 - y - y'l) + y'k-^c^''d,''c^d^{5 + k^jf; 

Dniy) = -r ' - y'7-'){5 - y'l) + y'k^c^d^c^''d^''{k-^6~' + ^^ f- 

B^iy) = y{6-' - y^j-^^k-^c^'^d^'^iS + k^j) + c^'' d^'' (k'^ 6-' + j-')}; 

Cniy) = -y{5 - yh){k'' c^^k-^ 5-^ + 7-^) + c^d^{5 + ^^7)} , 

where 5, 7 are defined by 

S .= do d^ =d^ d^ ,7 .= = -™ . 

Eliminating the y-factor in the equation of Bj^, we obtain an irreducible curve, which will still be 
denoted by Bj^ for the convenience of notations, but now with the equation: 

a(y2^-5)r?2 + by7? + c(yV'-5"') =0 , (55) 
where a, b, c are the parameters defined by 

a := k^c^d^{k-^5'^ + 7"i) + c^d^{5 + k^j) , 

b := k-^Co''do''c^d^i5 + fc^7)2 - k^c^d^c^''d^''ik-^5-^ + 7-^)2 , 
c := fc-^co ^^0^(5 + A:^7) + c^"" d^^ {k-^ + 

The curves Bf^ form a family of elliptic curves, depending on the four parameters, 5, 7, k^ Cq d^ , C2 d2 ■ 
And is a Z^-cover (branched) over i3^, with the covering transformation group containing t± in 

(g). For a generic /i, is a high-genus curve; indeed the genus is equal to 2N'^{N -1){N + 2) + 1. 

Now we are going to make a qualitative analysis on solutions of the Baxter's T-Q relation ([Tl| ) 
related to the discrete quantum pendulum (17) and SG model (p!8[). All the linear transformations 



4 

appeared in the expressions of (0) are operators of the vector space ® C^*. As D,C and UjS 
in the expression of Ti{x) are commuting operators, by { ^f^^ C)^ = !> the eigenvalue of TS{x) 



are still expressed in the form of (|2q ) with A depending on k; while for the discrete quantum 
pendulum and SG model, it becomes (p2[). By the expressions of D and C, it is not hard to see 

4 12 

that the common eigen-subspaces of q^* for the commuting operators D2^ ^ coc2 (j Jiave the 



C1C3 

N* 



dimension . The eigenspace decomposition of C * is denoted by 

(g) C^* = E„,„, , E„,„, C^' , 



i.rt'eZ 



JV 



where D^, -S^^C ^ act on E„ „/ by the multiplication of q"',q"'' respectively. By the relations of 

rC C0C2 ' 

?7j's and C^D in (p^), each E„^„' is stable under f/js. The operators UjS on E„^„' are determined 
only by those of Ui,U2, which form the Weyl algebra: U2U1 = UJU1U2 , Uf^ = U2 = 1- This 
implies the operator T| in (|l^ is determined by the representation E„.„' of the Weyl algebra 
on each sector, labelled by the eigenvalues of Tq^T^ corresponding to (n,n'). As the irreducible 
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representation of Weyl algebra is unique, given by the standard one on C'^, E„^„/ is isomorphic 
to the sum of A^-copies of the standard representation as Weyl algebra modules. In particular, the 
eigenvalues of — T| on the vector space En,n' are induced from the standard representation of the 
Weyl algebra ; each eigenvalue gives rise to eigenvectors in E„ „/. By ^ and (|5^, the Baxter 

4 ^ 

vacuum state \p) G C'^ for p € Cr is now defined by \p) = \po) \pi) \p2) Ips), where the 
vector \pj) in is given by the conditions: 



(Ob,- 



1 



{m\pj) 



-,2m-l 



+ xdj 



{m — l\pj) —S,j{6,j+ixcjq'^"^ — dj)k'^^Cjdj 



For a generic h, 
induces an isomorphism between ® 
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on the Baxter vacuum state, 



the evaluation of vectors of 

C*^ and a A^^-dimensional subspace of rational functions of 
C^. Through this isomorphism, E^^^/ gives rise a iV^-dimensional functional space of C^, denoted by 
e(E„ „/), with a Weyl-algebra-module structure induced from E„ „/. In the Baxter's T-Q equation 
([Tl|) on with A*(x) = Am,z(x) in (|^), the function Q{p) is the TJ(x)-eigenfunction in e(E„^„/) for 
(n, n') = (m, I), {N — m, N — I), with the multiplicity N. For the discrete quantum pendulum and 
SG model, by (|l^) and the A*(x) in Baxter's T-Q relation (1TT| ) is the reciprocal polynomial 



(32). To determine the eigen- functions Q{p) would require the understanding of its zeros and poles 
by using the expression of the Baxter vacuum state, which has been a difficult task at this moment. 
The possible role of elliptic function theory of Bj^ in the solutions of Baxter's T-Q relation on C^, 
and some further understanding on the eigenvalue problem of the models (|T^ (^8|) along this line, 
would be the core of our future work in this aspect. 



7 Conclusions and Perspectives 

We have studied the discrete quantum pendulum and discrete sine-Gordon model in the framework 
of quantum inverse scattering method. The diagonalization problem is governed by the Baxter's 
T-Q relation, which arises from the Baxter vacuum state on the spectral curve through a general 
scheme by using the transfer matrix for a fixed finite size L. We have demonstrated the role 



of algebraic geometry in the qualitative study of the Baxter's T-Q relation for L=3 in |13], and 
L=4 in this article. In both cases, they have shown an intimate relationship with integrable 
Hamiltonian spin-chains of physical interest. In this approach, one relies on the spectral curves, 
depending on the parameters encoded in the corresponding Hamiltonion expression. For generic 
parameters, the spectral curve, where the Baxter's T-Q relation is formulated, is a high-genus 
Riemann surface, as demonstrated in Sect. 6. However for both L = 3 and 4, the spectral curves 
possess a common feature that they form certain branched covers over elliptic curves. One might 
expect to employ the elliptic function theory to the solutions of Baxter's T-Q relation so that 
the algebraic geometry study could enrich our understanding of the corresponding Hamiltonian 
spectrum problem. Although this thinking is merely a speculation at present, a possible program 
along this line could be a challenging one, on which we hope to make progress in future. 

When the spectral curve degenerates into rational curves, where the geometry plays little role, 
we derive the polynomial formulation of Baxter's T-Q relation for a system of a finite size L in 
Sect. 3. We apply these results to the case L=4 in Sect. 4 for the setting of discrete quantum 
pendulum and sine-Gordon model. In these cases, an extra symmetry has naturally been imposed 
on the Baxter's T-Q polynomial equation; indeed it is governed by the reciprocal property of the 
equation. We present a detailed and rigorous mathematical derivation of the solutions in Theorem 
^. Surprisingly the conclusion on these polynomial solutions has been much in tune with the 
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one for L = 3 on the study of Hofstadter-type model (see Theorem 3 in |13|). Furthermore, 



the exact connection of the Baxter's T-Q polynomial equation with the Bethe ansatz technique in 
literature has been clarified in all these cases. The results obtained in this paper signal some further 
mathematical feature of the Baxter's T-Q polynomial equation, namely, the novel connection with 
certain (/-Sturm-Liouville problem at roots of unity = \. The facts discovered in this work 
could be served to demonstrate that a systematic mathematical theory embodied in the Baxter's 
T-Q polynomial equation (or algebraic Bethe Ansatz) would emerge in the study of g-difference 
operators. Accordingly, the relationship along this line is now under our consideration. 
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Appendix: Discrete Quantum Sine-Gordon Hamiltonian 

For the consistency of our notion of discrete quantum sine-Gordon Hamiltonian with the ones used 
in other literature, we make an identification of the discrete sine-Gordon integral in Q with the 
T/s of (|l|) in this paper. In Sect. 5 of @, the discrete quantum sine-Gordon Hamiltonian arises 
from the following commuting operators^ 

^(1)= U2Z2Ur^ + UiZiU^^ + V2h2hJVr\ 

^(2)= U2Z2U1Z1, 

where the lower index j = 1,2 indicates the site of operators in the same algebra generated by 
U,V, Z, subject to the relation UV = q 2 VU with Z the central element, and h is defined by 
h = /c^ + A;2q^U2Z. The sine-Gordon {SG) integral is defined by the operator 

One can show that 

V2h2htVi-i = q^V2Up2Vr^ + q^VaZr^UrX^i + WaUiZaZr^Ur^Vfi + k-^Vi^i^ ■ 
For the convenience in expressing H, we denote 

Wi := q5V2Zr^UrX-\ W2 := U1Z1U2 ^ , 

W3 :=q^V2U2Z2Vri , W4:=U2Z2Ur^ 

Note that 

W2W4 = . (56) 

We have 

fcV2UiZ2Z^iUf2Vf ^ = fcq^A(0)W3W^\ k-^y2yi^ = A:-iq^A(2)W4 ^Wi, 

''Here we use the sans serif type style, instead of the italic type style in for operators appeared in the right 
hand side of the expressions, for the purpose of less confusion on notations used in this paper. 
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hence 

^(1) = Wi + W2 + W3 + W4 + ytq's^ ^(o^WgW^i + k'^q"^ A^'^^^^^\Ni . 

With q2 = we identify the operators appeared in the above ^(•'^s with those in Tg^-s under the 
condition dQdid2d^ = 1, 

and then impose further constraints, 

1 ^ cMj^^^^ ^ ^(2) ^ ^ (57) 

A;c2 

By the equahties Vi = U3U2,Vi = U2U1 in (||), the 5G-integral # becomes -T^ in (|lD. Then 
(57) gives rise to the identification, 

W2 = W4 , = , 

or equivalently, 

_feC2_^i^_l ^ 03^1^2^^^ C1C3 ^ ^ 

cido(i3 ^ ^co ' k'^CQC2 

Note that the above relations are consistent with the relation U2U4 = C^^D in (^). Hence we 
obtain the relations (16) (18). 
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